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1. Introduction 



Perturbation theory is an indispensable calculational tool in particle physics. Meth- 
ods for perturbative calculations have been developed concurrenlty with the intro- 
duction of field theories for describing particle interactions. It is not surprising that 
we already have very efficient tools which confront experimental data at a quantitive 
level. 

It is evident, however, that the current methods are not suitable for comput- 
ing with sufficient accuracy all the required cross-sections at modern experiments. 
At the LHC or a future Linear Collider, for example, we must study a number of 
new complicated processes in the Standard Model or other theories. In addition, 
small effects arising at higher orders in perturbation theory will become significant 
in these experiments. It is important for such studies to improve or replace methods 
which require substantial human intervention. Ideally, we should develop automated 
methods applicable to every process, theory, and order in the perturbative expansion. 

Two types of computations are generally required for the evaluation of cross- 
sections and decay rates: loop integrations over the momenta of virtual particles, 
and phase-space integrations over the momenta of particles in the final state. At 
higher orders in perturbation theory both tasks are hard; this is primarily due to the 
large number of integrals which typically appear. Unfortunately, methods for the 
analytic computation of loop and phase-space integrals are complicated; it is usually 
unrealistic to attempt a brute-force computation for all terms in the matrix-elements. 
A solution to this problem is to construct algorithms which reduce the number of 
integrals to a few master integrals, and calculate directly the master integrals only. 

The method of integration by parts (IBP) for the reduction of loop integrals was 
introduced in [@, 0. Integrals which have common propagators (or, equivalently, 
belong to the same topology) satisfy linear algebraic identities. These identities can 
be derived with the IBP method and can be cleverly combined to produce reduc- 
tion identities to master integrals. Gehrmann and Remiddi introduced a new class 
of identities for scalar loop integrals due to their invariance under Lorentz transfor- 
mations Lorentz invariance (LI) identities are particularly useful for multiloop 
integrals with many external legs and massive propagators. Recently, the method 
of IBP and LI identities was extended to phase-space integrals that appear in the 
evaluation of total cross-sections 0, H and various differential distributions [P, 0, |]. 

Because of its conceptual simplicity, the IBP method was used to construct 
reduction algorithms for many classes of multi-loop integrals (see for example |TU 



Tl|, O, |13|). Nevertheless, the construction of such programs was laborious and a 



systematic approach to produce reduction identities for arbitrary topologies was not 
available; this was the main reason for the rather slow pace of multiloop calculations. 

This situation is now improved due to Laporta, who has proposed a fully au- 
tomated method for the reduction of generic loop amplitudes |T^. In contrast to 
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earlier approaches, his method does not attempt to derive reduction identities ap- 
phcaple to all the integrals of a topology. Instead, the aim is to reduce one-by-one 
the integrals by solving a large system of IBP/LI equations. This is achieved using 
Gauss elimination, after the IBP/LI system is ordered according to the complexity of 
the equations. Starting from the simplest one, each IBP /LI equation of the system is 
rearranged following a few algorithmic rules: the terms of the equation are assigned 
a relative weight for their complexity, and the most complicated term is then isolated 
on the left hand side. A recursive application of this procedure leads to expressions 
for complicated loop integrals in terms of master integrals. 

The algorithm proposed by Laporta has already been used in a variety of calcu- 
lations (for example in |l6l |l7l g g |g 0, |, 0, ||, |1 ) . However, we have 



found that its efficient implementation in a computer program is not trivial. The 
main difficulties arise from the fact that typical multiloop calculations, such as the 
ones mentioned earlier, require an enormous number of IBP/LI equations (10^ — 10^). 
In the process of Gauss elimination the algorithm can produce very large expressions; 
one must optimize for their efficient manipulation. 

In this paper, we provide a MAPLE 9 computer program (AIR) based on 
the method of [Q, for the Automatic Integral Reduction at higher orders in 
perturbation theory. The user should supply template IBP/LI equations for the 
integrals of a topology, optional information on the vanishing integrals of the topology 
and the master integrals (if known), and a small number of parameters controlling 
the treatment of large expressions. There is no need for advanced knowledge of the 
MAPLE platform. The input can be supplied with easy to modify text files, and 
AIR can be controled with very simple scripts. 

We believe that theorists who do not wish to invest in studying and implement- 
ing reduction methods, but need to study higher order effects in perturbation theory 
for various physical processes, will find this to be a valuable tool. We also hope that 
this publication will initiate some activity and exchange of ideas on practical issues 
concerning the implementation of reduction algorithms. In this program, we have im- 
plemented computational tricks wich were developed during practical computations; 
we hope our program will be improved from the experience of other users. 

The cost in computer resources grows rapidly with the complexity of the study 
process. It is inevitable that AIR will fail to solve arbitrarily large systems of equa- 
tions with large number of symbolic parameters (corresponding to kinematic scales, 
dimension, etc). However, we do expect AIR to be used for many applications in par- 
ticle phenomenology beyond the current state-of-the-art. For this purpose, we have 
included routines which minimize the number of computations during the reductions, 
mainly by keeping the number and the size of the actively processed expressions for 
Gauss-elimination to a minimum. 

In Section |^ we review the main features of the algorithm of Laporta using the 
massless one-loop box integral as a pedagogical example. In Section |^, we explain 
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the main features of AIR. In the rest of the paper we demonstrate the usage of 
AIR through examples. In Section | we reduce the massless one-loop box topology 
with no special algorithms activated for handling large expressions. In Section ^ 
we repeat the reduction by activating a "masking" algorithm which reduces the 
amount of computations during Gauss elimination by storing away the expressions 
which get reduced in terms of master integrals. In Section ^ we apply a different 
masking algorithm for very large integral coefficients. In Section |^ we show how 
to reduce topologies with a very large number of kinematic scales, by switching off 
simplification routines. We use the reduction of the massless one-loop pentagon 
topology as an explicit example. Finally, we present our conclusions in Section |. 



2. The reduction algorithm 

In this Section we present the reduction algorithm which is used in our program. A 
detailed description of the algorithm can also be found in Ref. |T^. Here we will 
present its main features using the massless one-loop box topology as an explicit 
example. 




We consider the class of integrals: 

_B(z/i, z/2, z/3, z/4) = f (fk y ^ I, (2.1) 

J [/c^J [(A; + pi)^J [(A; + pi2)^J ' [(fc + Pi23)^J 

where we have introduced the shorthand notation 'Pij,,,k = Pi + + ■ ■ ■ + Vk- The 
terms in the denominator are raised into positive or negative integer powers Vi. Zero 
powers correspond to scalar triangle and bubble integrals; negative powers correspond 
to triangle and bubble integrals with irreducible numerators. The external momenta 
are all taken to be light-like, v\ = J)^ = yi, = P123 = 0. These integrals arise in 
one-loop QCD amplitudes for 2 ^ 2 scattering processes (e.g. gg gg [p5[] ). 



It will be useful to know the values of the parameters for which the correspond- 
ing integrals vanish (tadpoles, scale-less bubbles). This information is not formally 
required; by solving the IBP equations one will eventually find that tadpoles, etc, 
are indeed vanishing. However, it is more efficient for the reduction to utilize the 
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fact that many terms in the IBP equations are zero. We find the following vanishing 
integrals: 



if 



or 



or 



or 



or 



e(z^i) + e(//2) + e(//3) + e(z/4) <2, 
e(//i) + e(//2) = 0, 
6(1^2) + e(z/3) = 0, 
6(1.3) + 6(1.4) = 0, 
e(//4) + e(//i) = 0, 



(2.2) 
(2.3) 
(2.4) 
(2.5) 
(2.6) 
(2.7) 



where we define Q{x) = 1 for x > and B(x) = for x < 0. 

Now we proceed to find algebraic equations for the integrals of the box topol- 
ogy. An easy way to derive such identities is the IBP method 0; we multiply 
the integrand with a loop or external momentum and differentiate it with the loop 
momentum. These total derivatives integrate to zero: 







d 



'dk, [kT [{k + Pi)T [{k + Pi2)T [{k + Pi23)T' 
where rj = k, k + pi, k + pi2, k + P123. We obtain four IBP identities: 



Ti : = [siyil+ + {d - 1^12334) - (^^il^ + ^^22+ + 1^44+ 
T2 : = \tu22+ + {d- //12344) - (1^1 1+ + /^22+ + z/33 

V22+ + Z/33+ + 



T. : 



sz.33^ + {d — i/i 



1234. J 



3 
4 
1 



B 
B 
B 
B 



(2.8) 

(2.9) 
(2.10) 
(2.11) 
(2.12) 



T4 : = [ti/44+ + {d- //12234) - (1^21+ + ;^33+ + Va4+) 2 
where the action of i+ (i^) increases (decreases) Vi by one in the integral S, e.g. 

3^5 = 5(i.i,I.2,Z/3±l,Z/4). (2.13) 

Products of operators have a straightforward interpretation, e.g. 

3+l-S = S(z/i-l,z/2,z/3 + l,z/4). (2.14) 

We have also used the shorthand: I'ijjk... = t^i + 2uj + Uk + . . ., and we define the 
usual Mandelstam variables s = Pi2> ^ = P23- 

The IBP Eqs. p^p.l2| and, optionally, the results of Eqs. p^p^ are sufficient to 
reduce any integral of the box topology to master integrals by using the algorithm 
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of Laporta fT^. In Ref. |jT^, the reader can find a detailed and complete description 
of the algorithm; here we intend to emphasize its salient features. The user is not 
required to have knowledge of the algorithm, however, some familiarity will be ben- 
eficial. We will describe the algorithm by tracing the first steps of our code when 
solving the box topology. For concretness, we will stop when the integral B(l, -1, 1, 
0) is reduced in terms of master integrals. 

• Seed generation: The program starts with the simplest list (z^i, 1^2, t's, ^^4) for 
which Blui, V2i T^i) is not vanishing. 

Si : {vuV2,Pz,Va) = (1,0,1,0). (2.15) 

will be our first 'seed' for generating identities from the topology IBP equa- 
tions Ti, . . . ,T4 (Eqs. p^p.l2|) , which we use as templates. 

• Identities generated from the template IBP equations: Our first identity is Eq. 
Ti substituting the values of {ui} found in Si: 

El : sS(2, 0,1,0) + (d- 3)5(1,0,1,0) = 0. (2.16) 

In El, we have already used our knowledge for the vanishing integrals of the 
topology (Eqs. M^). 

The above equation can be recast to express one of the two integrals in terms 
of the other. We would like to use such equations to express more complicated 
integrals in terms of simpler ones and, finally, in terms of the master integrals. 
It is therefore necessary to introduce criteria for the complexity of the integrals; 
the most complicated should receive first priority and will be isolated in the 
left hand side. 

• Integral priority criteria: We check on three parameters in order to isolate the 
most complicated integral. First we select the integrals with the largest number 
of propagators: 

Krop = T.^i''i)- (2-17) 

i 

If more that one integral has the maximum Mprop) we select the one with the 
largest sum of positive indices Vi 

N+ = Y.^{vi){v^-l). (2.18) 

i 

If more than one integral has the maximum values of Nprop and A^+, we select 
the one with the largest sum for the magnitutes of negative indices Pi 

M- = -Y.Q{-Vi)ui. (2.19) 

i 

If still there is an ambiguity, we randomly choose one of the integrals which 
has survived all three criteria. 
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• Rearranging the identities: Following the previous criteria, we find that B(2, 0, 
1, 0) is the most complicated integral in El. We then rearrange the identity to 
produce: 

El : B{2, 0, 1, 0) = -^—B{1, 0, 1, 0). (2.20) 

s 

We proceed, in the same manner, with the remaining identities T2,T3,T4 for 
Si. We obtain: 

E'l : 5(2,0,1,-1) = ((i-2)B(l,0,l,0)-5(l, 0,2,-1), (2.21) 

Q J 

Ef : S(l, 0, 2, 0) = -^B{1, 0, 1, 0), (2.22) 

s 

and 

Et : S(2, -1, 1, 0) = (d - 2)B(1, 0, 1, 0) - B{1, -1, 2, 0). (2.23) 

• Seed priority criteria: We have now processed all IBP equations for the first 
seed Si. It is therefore necessary to choose a new seed to obtain more identities. 
It is important to choose seeds that are most likely to produce equations coupled 
with the ones which have been processed earlier. For this purpose, we could 
select the seeds with the opposite priorities than the integral priorities, i.e. the 
seed with succecively minimum values for [A/'prop, A/"-, A/^]- However, the rules 
for choosing the seeds are mostly empirical and require some experimentation. 
In fact, the order for applying the criteria for minimum Af_ and jV^ can be 
judiciously chosen according to the class of integrals that the user needs to 
compute. For example, we could now pick cither (1,-1,1,0) or (2,0,1,0) as 
the next seed. Since our goal is to compute i?(l, — 1, 1, 0) it is better to choose: 

52:(i^i,i^2,i^3,i^4) = (1,-1, 1,0). (2.24) 

which in the IBP equations generates integrals with the same structure as 
the one we want to solve. Our program generates the seeds automatically; 
the user must provide the range of A/J,rop, A/"-, A/+ as input. It is relatively 
straightforward to decide the values for these parameters by inspecting the 
integrals that are required in the study process. 

• Substitutions and Gauss-Elimination: We now find a new feature in Eq. Ti for 
the seed 5*2: 

El : sB{2, -1,1,0) + {d- 2)B{1, -1, 1, 0) = 0. (2.25) 

The integral 5(2,-1,1,0) is isolated at the left hand side (Ihs) of a previ- 
ous equation {Ef). In such cases, we eliminate the known integral from the 
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equation. Substituting Eq. |2.23| , and applying the integral priority criteria, we 
have: 

El : B{1, -1,2, 0) = ^^-^B{1, -1, 1, 0) + {d - 2)5(1, 0, 1, 0) (2.26) 

• Back-substitution: E\ is solved in terms of an integral that can be substituded 
back to Ef. We can now see how the Laporta algorithm works in practice; 
by adding new equations to the already solved equations we form new sub- 
systems of coupled equations wich eliminate previously unknown integrals. In 
our example, performing the substitution of E^ into E\ we obtain, 

Et : S(2, -1, 1,0) = ^S(l, -1, 1, 0). (2.27) 

s 

We process two more equations for the seed S2: 

El : 5(2,-1,1,-1) = -t5(l, 0,1,0) + (rf- 1)5(1,-1,1,0) 

+5(1, 0,1, -1)-5(1, -1,2,-1), (2.28) 

and 

El : 5(1,-1,1,0) = -|5(1, 0,1,0). (2.29) 

In the last equation we have computed the integral that we wanted in terms of 
a simpler one: 5(1, 0, 1, 0). It is clear from the previous equations that B(l, 0, 
1, 0) is a master integral. 

In summary, the algorithm requires the succecive generation of identities with 
terms of increasing complexity. The newly added equations usually contain terms 
which are also found in equations generated at earlier stages; this produces small 
subsystems of coupled algebraic identities. A series of substitutions diagonilizes these 
algebraic subsystems and yields complicated integrals expressed in terms of master 
integrals. The algorithm is a clever implementation of Gauss elimination. It exploits 
the fact that Feynman integrals can be ordered according to very simple criteria. 

We demonstrated how the algorithm reduces a number of integrals belonging 
to the massless one-loop box topology. However, there was no step in the previous 
reduction that depended on the specifics of the topology. Therefore, this algorithm 
is suitable for the reduction of generic multiloop integrals or, more generally, of 
parametric functions which satisfy coupled algebraic identities (e.g. hypergeometric 
functions) . 

3. Features of AIR 

In this Section we describe the basic functions of our program. The program is in- 
cluded as a gzipped and tarred file in the source submission of the electronic preprint 
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for this paper, and can also be downloaded from Ref. [Q. It is convenient to unzip 
and untar the distribution file in a directory where AIR can be located permanently. 

/home> tar -zxvf air.tar.gz 
AIR/ 

AIR/main.map 
AIR/BOXA/ 

AIR/BOXA/input_boxa . map 
AIR/BOXA/script_boxa . map 
AIR/BOXB/ 

AIR/BOXB/input_boxb . map 
AIR/BOXB/script_boxb . map 
AIR/BOXC/ 

AIR/BOXC/input_boxc .map 
AIR/BOXC/script_boxc .map 
AIR/Pentagon/ 

AIR/Pentagon/ input .pentagon. map 
AIR/Pentagon/ script_pentagon . map 
AIR/Pentagon5/ 

AIR/Pentagon5/ input _pent agon . map 
AIR/Pentagon5/ script_pentagon . map 

The distribution includes the program file main. map, input files inpuL- ■ -.map, 
as well as MAPLE scripts script- ■ -.map for the example reductions in the rest of 
the paper. The program consists of MAPLE routines for generating seeds for the 
template IBP/LI identities, finding integral priorities, generating the IBP equations 
from the seeds, performing Gauss-elimination, masking large integral coefficients and 
reduced expressions, performing nested substitutions, and collecting the results. The 
function of the more important routines will be detailed in the following Sections. 

In typical multiloop computations, a large number of identites should be pro- 
cessed. A database system is therefore required to access, modify, and store the 
equations. We have implemented a rather simple database system, where each IBP 
equation is stored in a single file; the name of the file is determined by the integral 
on the Ihs of the equation. We also create separate auxiliary files which serve to 
point to the equations in the IBP system where a particular integral can be found. 
Our database system is very robust; however, it creates a rather extended tree of 
directories wich usually contain very short ASCII files. 

The program can perform very complicated multiloop reductions. It is often 
possible to simplify all the terms in the IBP equations as they get substituted and 
rearranged for Gauss elimination. However, if the topology depends on many kine- 
matic scales, or the IBP equations are loosely coupled (creating large subsystems of 
equations before they get diagonahzed) , or the values for Aiprop, A/± are large, it may 
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not be feasible to perform all simplifications within acceptable times or the available 
memory. We have implemented two algorithms to perform the reductions efficiently 
and reduce the amount of computations; the algorithms can be used independently 
or in conjuction. 

The first algorithm masks subexpressions which are reduced in terms of master 
integrals. The program detects the reduced expressions, stores them in files, and re- 
places them by an indexed symbol. Thus, the masked expressions are protected from 
subsequent manipulations during Gauss elimination. This feature is implemented 
recursively; whenever a new expression is written in terms of masked expressions 
and/or master integrals, the new expression is also masked. At the end of the re- 
duction, a scries of nested substitutions is required in order to rewrite the masked 
expressions in terms of the master integrals. We will discuss the required nested 
substitutions later; for now, we should note that the masking algorithm reduces sig- 
nificantly the amount of computations during the process of Gauss elimination. The 
algorithm requires that the master integrals of the topology are known. To deter- 
mine them, one can perform a less involved reduction for relatively small values of 
M± without using the masking algorithm. When the master integrals are found, 
the user can repeat the reduction for larger values of A/±, activating the masking 
algorithm. 

The second algorithm aims to reduce the size of the equations by masking all 
integral coefficients which are lengthier than a user-defined maximum value. During 
Gauss-elimination, however, some integral coefficients vanish; the elimination cannot 
take place if the coefficients contain masked expressions. To solve this problem we 
check numerically for cancelations. The user needs to provide as input, numerical 
values for all the parameters (kinematic scales, dimension) which enter in the sym- 
bolic expressions for the integral coefficients of the IBP/LI equations. The masking 
algorithm substitutes these numerical values and stores both the numerical result 
and the symbolic expression for the coefficients. The program determines if a coeffi- 
cient is zero by inspecting the numerical result, thus, avoiding complicated symbolic 
manipulations. To ensure that cancelations are not accidental, the program can per- 
form the numerical testing of the values of the coefficients for more than one choices 
of numerical values for the kinematic parameters and the dimension. The analytical 
value of the lengthy coefficients is computed at the end of the reduction, and only 
for the integrals that arc required for practical purposes. 

The purpose of the two algortihms is to remove complications from the symbolic 
manipulation of very large expressions. However, after Gauss-elimination we must 
still perform computations which were defered by using the masking algorithms, i.e. 
we must perform a series of nested symbolic substitutions for the used alias symbols 
in order to compute the masked expressions explicitly in terms of the kinematic 
parameters, the dimension, and the master integrals. It is possible to imagine that 
this additional computation is as difficult as using the program without the masking 
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algorithms, where all substitutions take place explicitly during Gauss elimination. 
This is not the case; usually, only a fraction of the total number of masked expressions 
is required for the integrals that appear in the matrix-elements of a physical process. 
For example, the integral 5(2,-1,1,0) of the previous Section does not appear in 
the computation of e.g. gg gg, however, it appears in the IBP equations. By 
using the masking algorithms, we avoid computing the masked expressions for many 
such integrals. 

The remaining nested substitutions can still be challenging for very complicated 
problems. One can resort to tricks such as expanding in the dimension parameter 
or in kinematic parameters (e.g. electron mass in Bhabba scattering), if this is 
justified from the physics of the process. However, this is rarely needed; there are 
many processes where we can perform the substitutions without giving up on a 
valid evaluation of the integral coefficients for all values of the kinematic parameters 
and the dimension. AIR includes general purpose routines for a straightforward 
computation of recursive substitutions; these routines attempt a brute-force symbolic 
simplification of all the intermediate expressions. It also provides the option to switch 
off simplification of expressions that exceed a maximum length or, if necessary, to 



transfer the most complicated substitutions to another platform, e.g. FORM |27 



In the following Sections we perform four example reductions which can serve as 
a tutorial for using AIR and its main features. A technical description of the AIR 
routines can be found in [|26|. The source code of the program is openly distributed; 
the users are free to modify it. The authors will be greatful to receive suggestions 
and constructive feedback. 



4. Reduction with no masking 

We now perform the reduction of the massless one-loop box topology. In this Section 
we do not activate any of the two masking algorithms. The input file and the 
corresponding script for the reduction can be found in the directory: 

/home/AIR/BOXA 

The input file for the reduction is named input_boxa.map. It contains variables 
which are used globally by AIR. These are: 

ibp_equations : = [ 

-nu3*B (nul-1 , nu2 , nu3+l , nu4) -nu4*B (nul- 1 , nu2 , nu3 , nu4+l ) 

-nu2*B (nul-1 , nu2+l , nuS , nu4) +nu3*s*B (nul , nu2 , nu3+l , nu4) 

+ (-nu3-nu2-2*nul+d-nu4) *B (nul , nu2 , nu3 , nu4) , 

nu4*t*B (nul , nu2 , nu3 , nu4+l) -nu3*B (nul , nu2-l , nu3+l , nu4) 

+ (d-2*nu2-nu3-nul-nu4) *B (nul , nu2 , nu3 , nu4) 

-nu4* B (nul,nu2-l,nu3,nu4+l)-nul*B(nul+l,nu2-l,nu3, nu4 ) , 
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-nul*B (nul+1 , nu2 , nu3-l , nu4) + (-nu2+d-nu4-nul-2*nu3) 

*B (nul , nu2 , nu3 , nu4) -nu4*B (nul , nu2 , nu3- 1 , nu4+l ) 

-nu2*B (nul , nu2+l , nu3-l , nu4) +nul*s*B (nul+1 , nu2 , nu3 , nu4) , 

-nul*B(nul+l,nu2,nu3,nu4-l)-nu2*B(nul,nu2+l,nu3, nu4- 1 ) 

-nu3*B (nul , nu2 , nu3+l , nu4-l) +nu2*t*B (nul , nu2+l , nu3 , nu4) 

+ (-nu3+d-nu2-nul-2*nu4) *B (nul , nu2 , nu3 , nu4) 

]: 



This is a list of template IBP identities (Eqs. p.9[ - |2.12| ) for the box topology. 
The program reads off some additional implicit definitions from the structure 
of the IBP equations. For example, it is now defined that the name of the 
topology is ' 'B' ' and the powers of the propagators are defined through the 
variables nul , nu2, nu3, nu4. 



ZERO_TOPOLOGIES:=[ 

ThetaF(nul) + ThetaF(nu2) + ThetaF(nu3) +ThetaF(nu4) < 2, 

ThetaF(nul) +ThetaF(nu2) =0, 

ThetaF(nu2) +ThetaF(nu3) =0, 

ThetaF(nu3) +ThetaF(nu4) =0, 

ThetaF(nu4) +ThetaF(nul) =0, 

NULL] : 



• This is a list of statements (Eq. p.3[ - p.7|) which undergo boolean evaluation when 
the propagator powers are substituted by integers. If any of the statements is 
true, then the corresponding integral is set to zero. 



MASTERS : = [] : 



• This variable activates the algorithm for masking reduced expressions to master 
integrals. It contains a list of known master integrals. In this example, we do 
not want to activate the algorithm; therefore we define the above variable to 
be an empty list. 



check_values : = [] : 



• This variable activates the masking algorithm for integral coefficients which ex- 
ceed a maximum value. It should contain numerical values for all the kinematic 
parameters and the dimension in the IBP/LI equations. For this example, we 
do not want to activate the masking algorithm and we set the variable to an 
empty list. 



MAXLENGTH:=1000: 
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• This variable is used by the masking algorithm for large integral coefficients. 
It defines the maximum length for a coefficient in order not to get masked. 
The length of an expression is measured by the number of characters in the 
expression and is determined by a MAPLE routine. The value of the variable 
is irrelevant if the check_values list is empty. 

MAXSIMPLIFY:=10~10: 

• This variable is used by the routines which perform nested substitutions for 
the masked expressions and the routines which display the final results. It 
sets a maximum length for the expressions that MAPLE is allowed to simplify. 
Larger expressions get substituted but not simplified. 

MAPLEMAXSUB:=10"10: 

• This variable is used by the routines which perform nested substitutions of 
masked expressions, and the routines wich display the results. It sets a maxi- 
mum length for the expressions that MAPLE is allowed to substitute. Larger 
expressions are not written explicitly and are kept masked. In this reduction 
we want all coefficients to be explicit, and we set the value of the variable to a 
practically unreachable value. 

VERBOSE :=FALSE: 

• This variable is used to display information about the progress of the program. 

If set to TRUE, the program outputs on the screen the seed that is processing, 
or, after Gauss-elimination is comleted, the index of the masked expression 
that is evaluating. 

The file script _boxa. map contains all the calls to AIR for reducing the box 
topology. We can run the script from the shell command line: 

/home/AIR/BOXA> maple script_boxa.map 

It will be more instructive for this first application to call the routines interactively 
from within the MAPLE platform. We fisrt launch MAPLE, 

/home/AIR/BOXA> maple 

and load the input for the topology and the main program: 

> currentdir ( ' ' /home/AIR/BOXA " ) : 

> read input _boxa. map' ' : 

> read ' '/home/AIR/main.map' ' : 
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The user should now perform the following tasks: 

• Seed generation: We create a list with sets of integers for deriving IBP equa- 
tions from the template equations of ibp -equations. This is accomplished by 
calling the routine 

SEEDGEN( "filename", maxtop, [minNprop,maxNprop], [minNmi- 
nus, mcLxNminus] , [minNplus, mcixNplus]) ; 

The first argument is the filename where the seeds will be written. The next 
argument, mcixtop, is a list with integers {i\ denoting the propagators raised 
to positive powers in the seed with the highest priority. For example, if 
we are inerested in reducing all the integrals of the box topology and all its 
subtopologies we should set maxtop = [1,2, 3, 4], indicating that all z/i, P2-, 1^3, 
can appear with positive values. If wc only required integrals of e.g. a t-channel 
triangle subtopology, we could set maxtop = [1, 3, 4], indicating that seeds with 
positive V2 do not need to be included in the reduction. The next three en- 
tries determine the range of Nprop-, -A/"-, and repsectively. These values are 
mostly empirical. A rule of thumb is that one should generate seeds wich in- 
clude the indices of the most comhcated integrals to be reduced, as well as the 
complete tower of integrals with lower priorities. To give a concrete example, 
we will generate seeds for the integrals that appear in the one-loop gg gg 
amplitude. Wc find integrals from all subtopologies (bubble and triangles) of 
the box topology, therefore we set 2 < Mprop < 4. We also find that all inte- 
grals have A/V = (there are no squared propagators in the amplitude), and 
can have up to 4 powers of irreducible numerators: < H- < 4. To emphasize 
differences in the running times for various algorithms of the program we will 
extend the latter interval to < jV- < 10. 

> SEEDGEN ("seeds. map " , [1, 2, 3, 4], [2, 4], [0, 10], [0, 0]); 

The routine produces a list of ordered seeds in the file /home/AIR/BOXA/ seeds .map. 

• Gauss elimination: We generate the IBP equations from the seeds and perform 
Gauss elimination by calling the routine 

Reducer("seeds_file", "monitor_file" , "RESULTS_DIR"); 

The first argument is the name of the file with the seeds for the reduction 
(as it was produced with SEEDGEN). The next argument is the name of a file 
which the program updates with the processed seeds; it serves to monitor the 
progress of the program. The last argument provides a directory path where the 
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program can deposit the database with the IBP equations. For our example, 
we type 



> Reducer C' seeds .map' ' , ' 'calc.map' ' , ''.''); 

It is useful to inspect the /home/AIR/BOXA directory. The program has cre- 
ated a number of subdirectories (B_l_3, B_l_3_4, etc) which correspond to the 
non-vanishing subtopologies of the box topology. In these subdirectories, AIR 
stores the IBP equations after they have been rearranged to isolate the most 
complicated integral at the Ihs of the equation. For example, Eq. |2.29| is stored 
in 

/home/AIR/B0XA/B_l_3/B_l_-l_l_0.map 

The first part of the file path /home/AIR/BOXA/ corresponds to the directory in 
the third argument of the Reducer command. The second part of the file path 
B_l_3 is created from the integers {i} for which the powers {ui} are positive. 
The ending of the file path is created from the indices {z/j} of the integral. 

Collecting the results: Essentially the reduction is now complete; the reduced 
integrals can be found in the files of the database tree for the IBP equations. 
Inspecting some of the equations we find that many integrals are reduced in 
terms of three master integrals: B(l, 0, 1, 0), B(0, 1, 0, 1), and B(l, 1, 1, 1). 
However, we also find integrals which are not fully reduced, e.g. the integral 
B(2, 1, 1, -11). It is usually observed that integrals with the same indices as in 
the seeds are fully reduced. Motivated from this observation, we have written 
a routine for collecting in a separate directory the seed integrals: 

tidy_list("seeds_file", "RESULTS_DIR" ); 

The first argument is a file with seeds (as generated by SEEDGEN). The last 
argument is the directory where the program has placed the results (as in the 
third argument for Reducer). For our example, we type: 

> tidy_list( "seeds. map " , ''.''); 

The routine has created a subdirectory, named RESULTS, wich contains the ex- 
pressions for the seed integrals only. For example, one can find the integral 
B(l, -4, 1, 1) in the file /home/AIR/B0XA/RESULTS/B_l_3_4/B_l_-4_l_l .map. 
We should note that when the masking algorithms are activated, the tidy_list 
routine also performs the required nested substitutions for expressing the co- 
efficients of the master integrals in terms of the kinematic parameters and the 
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dimension. However, it does so only for coefficients with smaller size than the 
MAPLEMAXSUB value. It is possible to achieve a fast execution of the tid.y_list 
routine if we choose a low value (100-500) of MAPLEMAXSUB. The task of comput- 
ing large coefficients is postponed further, and is performed by a new routine 
which we shall describe shortly. This routine is customized to disentagle the 
newest masking, which has a simpler structure than the masking of Gauss- 
elimination, very efficiently. 

• Reading the results interactively: A useful routine for reading the reduced in- 
tegrals from within the MAPLE environemnt is 

showJnt (integral, RESULTS_DIR); 

The first argument is the required integral, and the second is the directory 
of the reduction (as in Reducer). For example, 

> show_int(B(l, -4, 1, 1), 

outputs the expression for the required integral in terms of master integrals. 
This routine can be used for purposes of interfacing the results of the reduc- 
tion to other programs that users develop for calculating matrix-elements. We 
should note that if the masking of master integrals is activated and the value 
of MAPLEMAXSUB is small, show_int will return the wanted integral as a lin- 
ear combination of master integrals but with masked coefficients. The routine 
show_full_int displays the unmasked result. 

We have described the basic variables and routines of AIR by performing the reduc- 
tion of the box topology. It is worth noting that the total running time for the three 
main routines (SEEDGEN, Reducer, tidyJist) is approximately two minutes on a 
1.6GHz processor. 

5. Masking reduced expressions 

In this section we repeat the reduction of the one-loop massless box topology pro- 
viding the known master integrals as input to the program. The algorithm uses this 
information to find parts of expressions which are reduced to master integrals, and 
masks them. The reduction proceeds faster, having replaced the masked expressions 
by indexed symbols K{i). When the step of Gauss elimination is completed, only 
the indexed symbols in the expressions of the seed integrals need to be computed 
explicitly. 

The reduction is performed in the directory /home/AIR/BOXB. The program al- 
lows, in principle, activation of the masking algorithm for master integrals with- 
out changing directory. However, we recommend performing the reduction in a 
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new directory if new values for the variables ZERO_TOPOLOGIES, ibp_equations, 
MASTERS, check_values, and MAXLENGTH are required. This will prevent possible 
integral misidentifications caused from creating the database of IBP equations with 
inconsistent values for these variables. 

The input file input _boxb. map is modified to activate the masking algorithm for 
master integrals. In particular, the variable for the master integrals is not empty; we 
have defined 

MASTERS :=[ 
B(l, 0,1,0), 
B(0, 1,0,1), 
B(l, 1,1,1)] : 

These are the master integrals that we found in the first reduction. 
We now execute the MAPLE script in script_boxb.map: 

currentdir("/home/AIR/BOXB") : 
read " input_boxb . map " : 
read "/home/AIR/main.map" : 

SEEDGENC seeds. map", [1, 2, 3, 4], [2, 4], [0, 10], [0, 0]): 
Reducer ("seeds .map" , "calc.map", "."): 
tidy_list ("seeds .map" , "."): 

The script is executed in less than a minute; this is approximately 50% faster than 
without the masking algorithm. The routine Reducer has created a tree of subdirec- 
tories ICED/ICED#/KEXPR, where the masked expressions arc stored in files named 
kexpr_i.map. The expressions are replaced in the reduction by the symbol K{i), 
where i is an integer index. 

The masked expressions are defined recursively; it is therefore necessary to per- 
form a series of nested substitutions before we obtain their explicit form in terms 
of master integrals. This task is performed by the tidy_list routine, which stores 
explicit results for the masked expressions in the directory /home/AIR/BOXB/KMELT. 
The user is not required to know the details about the file structure for the masked 
expressions; the integrals in the directory RESULTS are fully evaluated in terms of 
master integrals and can be accessed as before. 

6. ]V[asking large integral coefficients 

In this Section we describe the function of the algorithm which masks large inte- 
gral coefficients in the IBP equations. We perform the reduction in the directory 
/home/AIR/BOXC, which contains the input file input_boxc.map and the appropriate 
MAPLE script script_boxc .map. We have modified the input variables: 
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MASTERS : = [] : 
check_values : = [ 
[s=l,t=-0.12,d=3.3] , 
[s=-1.2, t=12.2, d=42] , 
[s=1.82, t=-0.345, d=-28.1] 
]: 

MAXLENGTH:=50: 

In this run we have deactivated the algorithm for masking reduced expressions. In- 
stead, we have provided a hst of numerical values for the kinematic parameters and 
the dimension which activates the algorithm for masking the integral coefficients with 
length bigger that the value of the MAXLENGTH variable. The MAPLE commands for 
the reduction are collected in the file script_boxc .map, 

restart ; 

currentdir("/home/AIR/BOXC") : 
read " input _boxc .map" : 
read "/home/AIR/main.map" : 

SEEDGENC seeds. map", [1, 2, 3, 4], [2, 4], [0, 10], [0, 0]): 
Reducer ("seeds .map" , "calc.map", "."): 
tidy_list ("seeds .map" , "."): 

The script is executed in about 90 seconds; this is slower than previously. In general, 
masking large coefficients is not as fast as masking the reduced expressions to master 
integrals. However, in reductions of complicated IBP systems, such as in mulitloop 
crossed topologies, we have found that this algorithm is indispensable. In extremely 
complicated problems it is required that both masking algorithms are activated. 

In the directory ICED we find two new subdirectories: EXPR and NUM. The rou- 
tine Reducer saves the expressions for the masked large integral coefficients in the 
first directory. In the second directory, it stores the numerical values of the same 
coefficients for the choices of the parameters in check_values. Expressions in the 
directory EXPR are defined recursively through other masked expressions. However, 
in the directory NUM they are explicitly evaluated for the special input values of the 
parameters. The program sets an expression to zero in the IBP equations if its 
numerical values is zero for all input choices in check_values. It is important to 
provide "sufficiently random" lists of numerical values for the parameters in order 
to reduce the risk for accidental cancelations. It is important to note that floating 
point numbers in check_values are automatically converted to fractional numbers 
when used in the program. Therefore, the numerical evaluation of coefficients is 
exact (using integer arithmetics), avoiding complications due to rounding. 

Finally, the results are collected in the directory RESULTS using the routines 
tidy_list and show_int. We should stress, that the results 
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7. Reduction of one-loop pentagons 



In this section we perform the reduction of the one-loop pentagon with massless 
propagators. This is an example of a topology with many kinematic parameters. As 
is common in such topologies, the expressions for the reduced integrals are large. We 
will therefore use this example to demonstrate the options in AIR for dealing with 
large expressions. 




We consider the class of integrals: 

where P12345 = and pi = pi = pi = Pa = p\ = ^- The IBP equations for the 
pentagon topology are: 
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where we denote the invariant masses Sij = pf^. We also note that the topology 
vanishes if any three adjacent propagators are raised to non-positive powers. 

After a preliminary run with no masking algorithms we find 11 master integrals: 

• the pentagon integral V{1, 1, 1, 1, 1), 



• the box integals 1/(1,1,1,1,0), 1/(1,1,1,0,1), 1/(1,1,0,1,1), 1/(1,0,1,1,1), 
F(0,1,1,1,1), and 

• the bubble integrals 1/(1, 0, 1, 0, 0), 1/(1, 0, 0, 1, 0), 1/(0, 1, 0, 1, 0), 1/(0, 1, 0, 0, 1), 
and y(0,0,l,0,l). 
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Remarkably, by using the algorithm for masking reduced expressions, the step of 
Gauss elimination is very fast and can be performed for far more complicated inte- 
grals than the ones that are practically needed; it only takes a few minutes before a 
sufficiently large number of integrals for QCD five-point amplitudes is solved. How- 
ever, keeping the variable MAPLEMAXSUB in unreachable values (i.e. switching off the 
special routines for large coefficients, as we did in all previous examples), forces the 
routine tidy_int to carry the load of all substitutions and simplifications. The pres- 
ence of five independent Mandelstam variables and the dimension parameter makes 
the step of nested substitutions very hard for the integrals with highest priority. The 
routine can get stalled in computers with a memory smaller than 1GB. 

We set the variable MAPLEMAXSUB to a value of 100 in the input_pentagon.map. 
The reduction routine Reducer is not affected by the new setting, however, the nested 
substitution routine tidy_list runs differently; it does not perform simplifications on 
any expression with more than 100 characters. Instead, it replaces these coefficients 
with indexed aliases, {/[«]}• The routine does not process large coefficients, and 
finishes very quickly. The difficult substitutions for the masked expressions f[i] are 
performed with a new routine melt_all_f ("."). This routine is designed to work 
with the smallest possible memory consumption. Finally, the user can display the 
seed integrals explicitly written in terms of master integrals, the dimension, and the 
kinematic parameters, using the show_f ull_int routine. 

It is worth demonstrating some important technical details for the reduction. 
We perform the reduction in the directory /home/AIR/Pentagon, where we have 
placed the input file input_pentagon.map, and a script for executing the AIR rou- 
tines script_pentagon.map. As usual, in the input file we have provided the IBP 
equations, conditions for vanishing integrals, and the list of master integrals. We 
also set values for the variables, 

MAXLENGTH:=100: 
MAXSIMPLIFY:=1500: 
MAPLEMAXSUB: =100: 
VERBOSE :=FALSE: 

We now perform the reduction of integrals of the pentagon topology with jVL < 4, 
and J\f+ = 0. 

\home\AIR\Pentagon> maple script_pentagon.map 

The script is executed in approximately 10 minutes using approximately 40MB of 
memory. At this point it is worth making some observations about the results of 
the reduction. We launch a MAPLE window, load the AIR files for the pentagon 
topology, and read the result for one of the reduced integrals: 

> currentdir ( ' ' /home/ AIR/Pent agon ' ' ) ; 
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> read input .pentagon. map' ' ; 

> read ' '/home/AIR/main.map' ' ; 

> show_int(V(l, 1, 1, 1, -4), 

The result now contains masked expressions, e.g. f [67]. Their value is stored in the 
ICED directory tree. We can see the value of the masked expressions using the 
routine: 

show_f(i, RESULTS_DIR); 

which displays the variable f[i] for the integer i. We must also provide the reduction 
directory RESULTS_DIR For example, the value of f[12] can be retrieved by typing: 

> show.f (12, " ."); 
Moreover, the command 

> show_full_int(V(l, 1,1, 1,-4) , 

will return the full expression of the wanted integral. The commands show_int, 
show_full_int and show_f are very convenient to collect the results and export 
them to other platforms for further calculations. 

We proceed, next, with further reductions of integrals of the pentagon topology 
with M- < 5, and A/+ = 0. To perform this reduction we have to prepare a run using 
the same input file and requiring additional seeds generated in the script file, 

> SEEDGENC seeds. map", [1,2,3,4,5] , [2,5] , [0,5] , [0,0]); 

The script runs now in approximately 25 minutes making use of 65MB of memory. 
The coefficients involved are now larger, but they were computed very efficiently. 

We have performed a number of reductions for one loop hexagon and heptagon 
topologies and various two-loop topologies. For example, all loop and phase-space 
integral topologies in [|, ^ are fully reduced in less than 6 hours. The reduction of 
the massless double box topology with A/1 < 4, and A/+ = can be performed in 
about three days; the Reducer routine was running for approximately 7 hours, the 
tidy_list routine went through in one hour and ten minutes, while the melt_all_f 
routine worked for a couple of days. The cross-box massless box topology is reduced 
in approximately four days. In massive two-loop topologies, we have reduced integral 
topologies for the production of heavy quarks. The double-box topology with two 
massive external legs and a massive propagator (all carrying the same mass) is re- 
duced using the masking algorithms in about 32 hours while the nested substitutions 
are completed in about 20 days for A/1 < 4, and J\f+ = 0. 
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8. Conclusions 



We have presented a MAPLE program for Automated Integral Reductions in per- 
turbative calculations. Our program is based on an algorithm introduced by La- 
porta |T^, and uses the method of Gauss elimination for solving large systems of 



equations. The program can reduce generic loop or phase-space integrals or other 
functions (like hypergeometric functions) which satisfy coupled algebraic identities. 

The main obstacle in multiloop reductions is the large size of the symbolic ex- 
pressions. We have implemented two algorithms in our program which organize the 
reduction more efficiently and reduce the amount of computations. The routines 
mask reduced expressions or large integral coefficients. This enables solving the 
system of IBP equations without performing all substitutions explicitly. The compu- 
tationally intensive task of nested substitutions is performed only after the procedure 
of Gauss-elimination is completed, and for only a small fraction of expressions which 
appear in the final results. 

Reduction algorithms cannot be extended to arbitrarily large calculations due 
to finite computing resources. We believe, however, that many phenomenologically 
interesting problems are tractable using AIR. There are a few improvements that 
we expect to make in future releases, implementing a more fiexible database for 
storing the equations, and including more efficient algorithms for performing nested 
substitutions and simplifying large expressions. The methods described in 
appropriate for achieving this goal. 
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